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the  major  mode  supported  by  this  structure,  data  are  displayed  of  at- 
tenuation constant  and  phase  velocity  associated  with  a variety  of  design 
parameters.  For  a limited  number  of  cases  radiative  losses  and  reflec- 
tions at  circular  bends  in  the  line  are  shown,  as  are  suggested  feed 
and  termination  configurations.  Some  calculated  field  distributions 
are  also  presented,  as  are  some  estimates  of  the  reflection  expected 
from  an  idealized  (metallic  wire)  intruder.  Some  experimental  data 
taken  on  a scaled  model  are  Included  to  verify  certain  results  calculated 
from  theory. 

A second  type  of  line  investigated  theoretically  is  the  single 
coaxial  cable  whose  outer  shield  is  perforated  by  circumferential  rec- 
tangular slots,  each  one  narrow  in  longitudinal  extent,  subtending  an 
arc  of  arbitrary  length  and  spaced  uniformly  along  the  line's  extent. 

The  cable  has  an  imperfectly  conducting  circular  center  conductor  sep- 
arated by  a lossy  dielectric  from  the  outer  shield  which  is  assumed 
to  be  perfectly  conducting  and  infinitesimally  thin.  This  configuration 
is  enclosed  in  a circular  sheath  of  lossy  dielectrirami  the  whole  is 
assumed  to  be  suspended  above  or  buried  within  a lossy  Nat  earth.  For 
one  geometry  corresponding  to  a commercially  available  version  of  this 
so-called  "leaky"  coaxial  cable  (Andrew  Corporation's  Radiax  cable) 
data  are  presented  of  the  attenuation  constants  and  phase  velocities 
of  the  major  modes  found  to  exist  on  the  structure.  A few  examples 
of  the  corresponding  field  configurations  are  also  given. 

As  a prelude  to  work  on  the  leaky  coaxial  cable,  some  effort  was 
expended  to  develop  the  theory  of  the  Goubau  line  over  the  flat  earth 
as  well  as  that  of  the  coaxial  cable  with  360°  circumferential  slots 
over  the  flat  earth.  In  addition,  although  tangential  to  the  main 
thrust  of  the  contract,  theory  was  developed  for  the  leaky  coaxial  cable 
centered  in  a cylindrical  tunnel  through  earth.  This  latter  theory 
could  easily  be  changed  to  accommodate  a cable  located  eccentrically 
in  a tunnel. 

The  report  concludes  with  a discussion  of  the  utility  of  the  lines 
as  intruder  sensors  and  suggestions  for  future  work. 
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1.  This  contract  was  to  investigate  guided  wave  structures  of  about 
one  mile  long  which  might  be  used  as  intruder  sensors.  The  guided  wave 
transmission  line  was  to  operate  in  the  presence  of  a lossy  earth.  The 
predominant  structure  examined  was  the  Vagi  Transmission  line.  It 
consisted  of  uniformly  spaced  imperfectly  conducting  thin  rods  (dipoles) 
of  equal  length  held  above  and  perpendicular  to  a lossy  earth.  The 
effect  of  bends  were  included  in  the  study.  The  attenuation,  phase 
velocity  and  surface  wave  field  were  determined  as  a function  of  the 
various  parameters  of  the  line.  The  effect  of  a simulated  intruder 

was  also  determined. 

2.  It  was  found  that  a Yagi  line  was  not  sensitive  enough  t detect 
an  intruder,  although,  improvements  do  seem  possible. 

3.  Generally,  much  was  learned  of  the  Yagi  line  and  the  influence  of 
its  parameters.  Enough  knowledge  has  been  attained  so  as  to  be  able  to 
consider  making  use  of  the  line  in  conjunction  with  other  structures. 

The  treatment  here  is  quite  thorough  and  permits  speculation  for  future 
work. 

t k ...  ( c C 
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I.  INTRODUCTION 


In  the  not-distant  future  It  is  expected  that  critical  permanent 
areas,  such  as  nuclear  power  plants,  nuclear  waste  sites,  dams,  water 
supplies,  prisons,  etc.,  as  well  as  semi -permanent  areas  such  as  for- 
ward located  airstrips,  docking  sites,  etc.,  which  are  vulnerable  to 
noncooperative  intruders,  will  require  protection  by  one  or  several 
means,  from  physical  barriers  to  electronic  sensors.  With  this  poten- 
tial need  in  mind,  the  present  contract  was  undertaken  to  investigate 
and  evaluate  a few  electromagnetic  guided  wave  structures  which  appear 
to  be  suitable  for  use  as  intruder  sensors  for  barriers  of  length  up 
to  one  mile.  Of  necessity,  any  such  structure  must  operate  in  the  pres- 
ence of  lossy  earth  and  theory  of  the  resulting  attenuation  of  the  guided 
wave  must  take  this  source  of  loss  into  consideration.  Additionally, 
it  was  suggested  that  the  guided  field  structure  should  be  sensitive 
to  an  intruder  within  a longitudinally  oriented  region  approximately 
2 meters  wide  and  3 meters  high  and  be  such  that  it  cannot  be  penetrat- 
ed from  beneath  without  detection.  Frequency  of  operation  was  speci- 
fied to  lie  in  the  range  50-500  MHz,  particularly  below  100  MHz  where 
data  indicate  a broad  resonance  in  standing  human  scatterers.  This 
report  summarizes  the  results  of  theoretical  analyses  of  two  types  of 
guided  wave  structures  - the  Yagi  transmission  line  and  the  "leaky" 
coaxial  cable.  Peripheral  results  are  also  presented  for  the  Goubau 
line. 


The  Yagi  (or  Yagi-Uda)  transmission  line  consists  of  uniformly 
spaced,  imperfectly  conducting  thin  rods  of  identical  length  suspend- 
ed a uniform  height  above  and  perpendicular  to  a lossy  flat  earth*'7. 

For  the  major  mode  supported  by  this  structure,  data  are  presented  of 
attenuation  constant  and  phase  velocity  associated  with  a variety  of  design 
parameters.  For  a limited  number  of  representative  cases  radiation 
losses  and  reflections  at  90°  circular  bends  in  the  line  are  shown, 
as  are  suggested  feed  and  termination  configurations.  Some  calculated 
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field  distributions  are  also  presented,  as  are  estimates  of  the  reflec- 
tion expected  from  an  idealized  (metallic  wire)  "intruder".  Limited 
experimental  data  taken  on  a scaled  model  are  included  to  verify  cer- 
tain results  calculated  from  theory.  Improvements  in  the  theory  are 
also  suggested  and  discussed. 

18  o6 

The  leaky  coaxial  cable  is  being  investigated  by  several  groups 
as  a potential  intruder  sensor  and  their  experimental  results  indicate 
that  it  is  a promising  candidate.  At  the  start  of  the  present  contract 
we  or  the  sponsor  know  of  no  theoretical  model  developed  for  the  par- 
ticular type  of  cable  investigated  here  from  which  numerical  results 
could  be  calculated  for  propagation  constants  and  field  structure. 

It  was  considered  useful  that  such  a theoretical  model  be  available 
to  help  explain  some  of  the  observed  characteristics  of  leaky  coaxial 
cable  and,  should  it  be  incorporated  into  sensor  systems,  to  provide 
some  guidance  in  the  design  of  the  cable  itself  and  its  deployment  geo- 
metry. Thus,  some  effort  was  expended  in  the  latter  half  of  the  present 
contract  to  develop  a theory  for  the  single  leaky  coaxial  cable  in  the 
presence  of  a lossy  flat  earth.  Although  this  cable  is  usually  used 
in  a parallel  pair  configuration,  time  considerations  forced  us  to  leave 
this  interesting  extension  undone. 

The  geometry  considered  consists  of  a single  coaxial  cable  whose 
imperfectly  conducting  circular  center  conductor  is  separated  by  a lossy 
dielectric  from  the  outer  shield  which  is  assumed  to  be  perfectly  con- 
ducting and  infinitesimally  thin.  This  outer  shield,  rather  than  braided, 
is  perforated  by  circumferential  rectangular  slots,  each  are  narrow 
in  longitudinal  extent,  subtending  an  arc  of  arbitrary  length  and  spaced 
uniformly  along  the  line's  extent.  This  configuration  is  enclosed  in 
a circular  sheath  of  lossy  dielectric  and  the  whole  is  assumed  to  be 
suspended  above  or  burled  within  a lossy  flat  earth.  For  one  geometry 
corresponding  to  a commercially  available  version  of  this  line  (Andrew 
Corporation's  Radiax  line  (27))  data  are  presented  of  the  attenuation 
constants  and  phase  velocities  of  the  major  modes  found  to  exist  on 
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the  structure.  A few  examples  of  the  corresponding  field  configurations 

are  also  given.  As  a special  case  which  is  interesting  in  Its  own  right, 

the  leaky  coaxial  cable  was  assumed  to  reside  along  the  centerline  of 

a cylindrical  tunnel  through  lossy  earth.  If  necessary  the  theory  can 

readily  be  modified  to  accommodate  an  eccentric  cable  and  thus  form 

the  basis  for  designing  an  open  communication  system  in  tunnels  - a 

00.00 

subject  of  interest  in  recent  literature  . 

The  report  closes  with  a discussion  of  the  utility  of  the  lines 
considered  here  as  intruder  sensors  and  suggestions  for  future  work. 

II.  THE  YAGI  TRANSMISSION  LINE 

A.  Synopsis  of  Theoretical  Approaches 

Although  most  theoretical  work  under  this  contract  centered  on 
the  straight  Yagi  line  of  infinite  extent  because  this  geometry  pre- 
sented a tractable  boundary  value  problem  in  classical  format,  certain 
questions  posed  by  the  sponsor,  such  as  those  Involving  radiative  losses 
at  bends,  feed  and  termination  configurations,  and  reflection  from 
an  intruder,  could  best  be  answered  for  a finite  Yagi  line  using  a com- 
puter-oriented approach  such  as  that  provided  by  the  method  of  moments. 

In  this  section  we  summarize,  in  varying  detail,  the  two  approaches. 

1.  The  Infinite  Line 

In  this  section  expressions  are  developed  for  the  surface-wave 
fields  and  the  transcendental  equation  whose  root  determines  the  com- 
plex propagation  constants  of  the  surface  wave  supported  by  an  inf i- 

1 2 

nitely  long  Yagi  line  in  free  space  and  over  the  flat  earth.  * 

Figure  1 illustrates  a periodic  array  of  dipoles  in  free  space. 

Each  dipole  is  a straight  wire  with  length  l - 2h,  radius  "a"  and  con- 
ductivity a.  Each  wire  is  parallel  with  the  z axis,  and  s denotes  the 
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REGION  3 


Figure  1.  Periodic  array  of  dipoles  in  free  space. 

spacing  between  the  axes  of  neighboring  wires.  The  time  dependence 
eiwt  is  understood.  The  surface  impedance  Z$  is  defined  in  terms  of 
the  field  on  the  surface  of  dipole  number  0 as  follows: 

Ex(a,x)  * Zs  H^(a,z)  for  -h  < z < h . (1) 

Tmpere's  law  relates  to  the  current  distribution,  so  that  Equation 
(11  can  be  written  as  follows: 

Ez(a’x)  = Zs  (2) 

where  IQ(z)  represents  the  current  distribution  on  dipole  0.  A valid 
relation  can  be  obtained  by  multiplying  both  sides  of  Equation  (?)  by 
IQ(z)  and  integrating.  This  leads  to  the  following  relaxed  boundary 
condition: 

?h  >o<z>  Ez  Ez'z>  dz  ■ Ki  'o<z>  dz  (3) 
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(4) 


E 


z 


Ez  = Eo 


l + E.J 


n=l 


(5) 


where  Ep  denotes  the  free-space  field  of  dipole  n and  denotes  the 
field  reflected  from  the  flat  earth. 


Let  IQn  denote  the  current  at  the  center  of  dipole  n.  When  a sur- 
face wave  propagates  on  the  dipole  array,  the  dipole  currents  are  given 
by 


I = I enYS 
on  oo 


(6) 


where  y denotes  the  propagation  constant  for  the  surface  wave.  From 
Equations  (l)-(6),the  terminal  voltage  at  dipole  0 is 


Vfy)  - ZQ  + ZQ  + 2 ^ ^ cosh(nys)  = 0. 


m 


The  term  ZQ  vanishes  for  a dipole  array  in  free  space.  The  mutual  im- 
pedance between  dipoles  0 and  n in  free  space  is  defined  by 


Zn  = 


-1 


*00  *on 


;h  Enfz)  dz 


(8) 


If  the  dipole  length  l does  not  exceed  X/2,  a suitable  trial  func- 
tion for  the  current  distribution  on  dipole  n is 


- 'on 


sin  k(h  - )z  |) 
sin  kn 


(9) 


where  k = Corresponding  to  equations  (8)  and  (9),  the  mutual 

impedance  Zp  is  given  by  Schelkunoff  and  Frlis39  for  parallel  filamen- 
tary sinusoidal  dipoles  with  spacing  p = |ns  |.  The  self  impedance  of 
dipole  0 (in  free  space)  is 
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(in) 


[kh  - sin  (kh)  cos(kh)]  Z 

7 = / + * 

0 0L  2tt  ka  sinZ(kh) 

where  ZQ0  denotes  the  self  impedance  of  a perfectly  conducting  dipole 

in  free  space.  To  calculate  Z , we  consider  the  mutual  impedance  of 

oo  r 

filamentary  sinusoidal  dipoles  with  spacing  p = a.  The  term  ZQ  repre- 
sents the  contribution  to  V(y)  arising  via  reflection  from  the  flat 
earth,  and  it  is  considered  later. 

In  Equation  (7),  V(y)  denotes  the  voltage  appearing  across  the 
short-circuited  terminals  of  dipole  0 with  a terminal  current  IQ0  of 
one  ampere.  This  voltage  must  vanish,  and  Equation  (7)  is  the  tran- 
scendental equation  for  the  system.  The  propagation  constant  y for 
the  surface  wave  is  determined  by  searching  for  a root  of  V(y).  As 
usual,  we  let  y = a + ,jB.  For  a perfectly  conducting  array  over  a loss- 
less dielectric  half-space,  the  attenuation  constant  a vanishes  and 
the  series  in  Equation  (7)  converges  slowly.  When  the  dipoles  and/or 
the  earth  have  finite  conductivity,  a is  positive  and  the  series  di- 
verges. To  improve  this  situation,  we  add  and  subtract  (in  Equation 
(7))  the  following  asymptotic  form  of  the  mutual  impedance: 

Zn  ~ (D/n)  e"jkns 

o 

n _ j (1  - cos  kh) 
u 

7iks  sin  kh 

where  n = J uQ/^0-  This  leads  to  the  following  transcendental  equation: 

V(y)  = ZQ  + z£  - D 8,n[l  - exp(ys-jks)]  - D &n[l  - exp(-ys-jks)] 

+ 2 f [7n  - (D/n)e'jkns]  cosh(nys)  = 0.  (13) 

n=l 

When  a and  Z vanish.  Equation  (13)  reduces  precisely  to  the  transcen- 
dental  equation  used  earlier  for  a perfectly  conducting  dipole  array 
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(11) 

(12) 
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in  free  space.  For  this  and  other  reasons,  it  is  convenient  to  refer 

to  "analytic  continuation"  as  the  basis  for  Equation  (13).  When  a = 

0,  the  series  in  Equation  (131  converges  rapidly  and  the  following  number 

14 

of  terms  is  sufficient  : 

N = 1 + 17  h/s  . (14) 

When  a is  positive,  we  still  use  Equations  (13)  and  (14).  Newton-Raphson 
iteration  is  convenient  for  locating  the  root  of  Equation  (13). 

In  deriving  Equation  (13)  we  used  the  following: 

00  Qnz 

£V  = -*n(1-e)  (15) 

where  z denotes  a complex  constant.  Equation  (15)  can  be  found  in  the 
literature  for  the  special  cases  where  z is  imaginary  or  negative  real. 

In  the  general  case  it  can  be  derived  formally  by  integrating  (with 
respect  to  z)  both  sides  of  the  following  summation  formula: 

l enz  = ez/(l  - ez)  . (16) 

1 

If  the  hyperbolic  function  in  Equation  (7)  is  written  as  the  sum 
of  two  exponentials,  then  Equation  (7)  contains  the  sum  of  a divergent 
series  and  a convergent  series.  The  divergent  series  is  given  by 

vl(Y)  = | Zn  e"^  . (17) 

From  Equations  (11),  (15)  and  (17), 


V^y)  = - D An  [1  - e(Y‘jk)s]  + | [zn  - (D/n)e"jknsJ e"Ys.  (18) 

When  a is  positive,  the  series  eventually  diverge  in  Equations  (17) 
and  (18),  although  they  appear  to  be  converging  when  n is  small.  There- 
fore, results  obtained  with  the  transcendental  Equation  (13)  are  not 
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considered  reliable  until  they  are  confirmed  by  other  methods.  Confir- 
mation has  been  obtained  by  several  methods. 


In  method  1,  we  use  Equations  (11)  and  (17)  to  obtain 


V 


1 


nys 


+ 


e(y-jk)ns 

n 


(19) 


where  M is  large  enough  to  justify  the  asymptotic  form.  If  as  and 
( B-k ) s are  small,  the  last  series  in  Equation  (19)  is  nicely  approxi- 
mated by  an  Integral  as  follows: 


vx  = X z enYS  + D / 

1 l n m+% 


* a(y-jk)sx 


dx 


(20) 


This  integral  is  recognized  as  the  exponential  integral  E^(z)  where 
z = - (M+M (y-jk)s.  Of  course,  the  integral  in  Equation  (20)  diverges, 
and  we  may  replace  it  with  E^(z)  only  If  we  are  willing  to  accept  the 
implied  deformation  of  the  path  of  integration.  When  as  and  (B-k)s 
are  small,  the  results  obtained  with  method  1 (Equation  (20))  agree 
closely  with  those  obtained  with  analytic  continuation  (Equation  (13)). 


When  method  1 begins  to  fail,  method  2 becomes  effective.  In  method 
2 we  consider  a large  but  finite  linear  array  of  parallel  dipoles  with 
uniform  length  and  spacing.  We  insert  a one-volt  generator  at  the  center 
of  the  dipole  at  one  end  of  the  array,  invert  the  impedance  matrix  and 
thus  solve  for  the  currents  on  all  the  dipoles.  Since  the  impedance 
matrix  has  the  Toeplitz  form,  it  takes  only  a few  seconds  to  carry  out 
this  operation  for  an  array  of  200  dipoles.  Now  the  current  I is 
known  at  the  center  of  each  dipole.  On  the  dipoles  near  an  end  of  the 
array,  the  current  on  one  dipole  may  bear  little  relation  to  the  currents 
on  the  neighboring  dipoles.  For  a small  group  of  dipoles  at  the  center 
of  the  array,  however.  It  is  reasonable  to  assume  the  currents  are  re- 
lated as  follows: 


I 


ny-,  s -ny.s  ny,s  -ny,s 

Irtn  = ci  e + C,  e + C,  e ' + C-  e 6 
on  1 2 3 4 


(21) 


Thus  we  are  postulating  two  modes  propagating  in  each  direction  along 
the  array.  There  are  four  unknown  coefficients  Ci  and  two  unknown  prop- 
agation constants.  By  using  the  recursion  relations  for  the  exponen- 
tial and  hyperbolic  functions,  it  turns  out  that  we  can  easily  solve 
for  y^  and  y?  if  the  currents  IQn  are  known  on  any  group  of  six  con- 
secutuve  dipoles.  In  this  manner  we  calculated  y^  and  y^  for  an  array 
of  190  dipoles  with  various  lengths  and  spacings.  In  all  cases  where 
we  were  able  to  fund  a root  of  Equation  (13),  it  turned  out  that  y^ 
was  associated  with  the  surface  wave  of  interest  here.  Furthermore, 
when  as  and  (B-k)s  are  not  too  small,  method  2 yields  values  of  y^  show- 
ing close  agreement  with  those  obtained  with  analytic  continuation. 

The  values  of  y?  indicate  a wave  propagating  essentially  with  the  speed 
of  light  and  with  a comparatively  large  attenuation  constant.  Although 
this  wave  is  not  well  understood,  it  may  be  the  free-space  radiation 
emanating  from  each  end  of  the  finite  array.  Waves  of  this  type  would 
presumably  propagate  at  the  speed  of  light  and  would  not  possess  a true 
attenuation  constant.  Instead,  they  would  behave  like  (l/R)e~JkR  where 
R denotes  the  distance  from  a "radiation  center".  Over  the  range  of 

six  consecutive  dipoles  at  the  center  of  the  array,  however,  the  func- 

-yR 

tion  1/R  can  be  approximated  by  the  function  e . 


It  was  found  that  method  2 fails  if  the  last  two  terms  are  deleted 
in  Equation  (21).  This  indicates  that  two  distinct  modes  are  prominent 
even  in  the  central  region  of  an  array  of  190  dipoles. 


Finally,  conservation  of  energy  was  employed  to  check  the  atten- 
uation constant  obtained  via  Equation  (13).  The  check  proceeded  as 
follows.  We  considered  a surface  wave  propagating  in  the  negative  y 
direction  on  a periodic  array  of  dipoles  in  free  space  as  in  Figure 
1.  Using  the  field  expressions  presented  here,  we  Integrated  the  Poynting 
vector  over  the  plane  surface  y = y^  to  determine  the  time-average  power 
transmitted  in  the  f-y)  direction.  Conservation  of  energy  requires 
that 
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Pt(s/2)  - Pt(-s/2)  * Pd 


(22) 


where  Pt(y)  denotes  the  time-average  power  transmitted  and  Prf  denotes 
the  time-average  power  dissipated  in  dipole  0.  Using  the  value  of  8 
obtained  via  Equation  (13),  we  calculated  the  various  quantities  in 
Equation  (22)  and  used  Newton-Raphson  iteration  to  adjust  a until  Eq- 
uation (22)  was  satisfied.  It  was  found  that  the  value  of  a obtained 
in  this  manner  agreed  closely  with  that  obtained  via  Equation  (13). 

In  view  of  the  results  obtained  with  method  1,  method  2 and  con- 
servation of  energy,  we  concluded  that  the  method  of  analytic  continu- 
ation (Equation  (13))  Is  valid  even  though  its  derivation  is  question- 
able. 


As  is  done  in  Equation  (4),  It  Is  convenient  to  express  the  surface- 
wave  field  as  the  sum  of  the  Incident  field  (generated  by  the  dipole 
currents  radiating  in  unbounded  free  space)  and  the  reflected  field 
arising  via  reflection  from  the  air-earth  interface.  We  are  now  prepared 
to  develop  suitable  expressions  for  the  incident  field. 

It  must  be  kept  in  mind  that  we  require  incident  field  expressions 
in  a form  that  will  permit  application  of  the  boundary  conditions  at 
the  air-earth  interface.  The  Poisson  sum  formula  may  be  thought  to 
be  useful  for  this  purpose.  Difficulties  arise  however  when  one  at- 
tempts to  apply  the  Poisson  sum  formula  to  the  field  of  the  surface 
wave.  Although  the  field  of  one  sinusoidal  dipole  has  a simple  closed 
form,  a divergent  series  results  from  summing  over  all  the  dipoles  in 
the  array.  However,  Floquet's  theorem  suggests  an  alternative  approach. 

To  simplify  the  field  expressions,  let  us  replace  the  round  wires 
with  thin  metal  strips  lying  in  the  yz  plane.  Then  the  volumetric  cur- 
rent density  for  dipole  n is  given  by 
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Jn  = z eyns  6(x)  P(y-nx)  IQ(z)/w  amperes/m2 


where  y denotes  the  propagation  constant  for  the  surface  wave  on  the 
array  of  dipoles  over  the  flat  earth,  and  P(y-ns)  denotes  a pulse  func- 
tion with  unit  height  and  width  w,  centered  at  y = ns.  So  that  the 
field  of  the  strips  will  closely  resemble  the  field  of  the  round  wires, 
we  choose  w = 4a  for  the  width  of  the  strips.  If  w is  much  smaller 
than  s,  we  can  replace  eyns  with  eYy  in  Equation  (23)  because  the  two 
functions  are  essentially  equal  in  the  small  region  where  P(y-ns)  does 
not  vanish.  Using  the  Fourier  transform,  we  can  express  Equation  (23) 
as  follows: 


z IQ(z)  eyy 


sin(pw/2) 

P 


e-jpns  ejpy 


dp  / cos(fx)  df. 
o 


By  summing  on  n and  using  Equations  (3-41)  and  (3-44)  in  Reference  40, 
we  obtain  the  current  density  of  the  periodic  array  of  strip  dipoles: 

I (z)eyy  °°  oo 

J = r l S (y)  J cos(fx)df  (25) 


c / \ _ sin(imrw/s)  J2rmry/s 
nry'  rmrw/s 


where  IQ(z)  is  given  by  Equation  (9)  and  SQ(y)  = 1. 

Using  Equation  (25),  it  is  not  difficult  to  show  that  the  free- 
space  field  of  the  periodic  array  of  dipoles  is: 


= - A eYy  l Sm(y)  / f sin(fx)  Q(z)  df 
m=-oo  o 


(28) 


Ey  * A £ (^mn/s)  Sjy)  / cos(fx)  Q(z)  df 

r 

. J - A e^y  J Sm(y)  / (k2+g2)  P(z)  df,  in  Regions  1,3  (29 a) 

Ez 

2kA  e^y  J Sm(y)  / [e'glzlcos(kh)-e'gh  cosh(gz)]^s-f^  df, 

in  Region  ? (29b) 


A - uu 

" ks  sin  kh 
g2  - f2  - (Y+j2itm/s)2  - k2 

P(z)  = / sin  k(h  - |z'|)  e“glz_z  I dz' 

-h 

ll 

Q(z)  = / sgn(z-z')  sin  k(h  - |z'|)  e"gtz‘z  I dz'. 


(30) 


(31) 


(32) 


(33) 


In  regions  1 and  3 of  Figure  1,  where  |z|  > h, 

P(z)  3 2k  e”glzl  (cosh  gh  - cos  kh)/(k2  + g2)  (34) 

Q(z)  = P(z)  sgn(z)  . (35) 

In  region  2,  where  |z|  < h, 

P(z)  * 2 £g  sin  k(h-|z|)  + k cosh(g| z| )e"gh  - k cos(kh)e“9l zt] /(k2+g2) 

(36) 

Q(z)  * 2k[cos  k(h-|z|)  - s1nh(g|  z|  )e‘gh  - cos(kh)e_gl  zl] 

k +gZ  • 


The  magnetic  field  intensity  is  found  from  curl  T,  and  it  is  observed 
that  hJ  = 0.  It  has  been  verified  that  the  given  field  satisfies  the 
appropriate  wave  equation  in  each  region  and  continuity  of  the  tangential 
field  components  across  the  boundaries  between  adjacent  regions. 

Now  that  suitable  expressions  are  available  for  the  incident  field, 
we  are  in  a position  to  consider  the  fields  reflected  from  the  air-earth 
interface. 

Figure  2 illustrates  a periodic  array  of  vertical  dipoles  over 
the  flat  earth.  The  field  incident  on  the  air-earth  interface  is  given 
by  Equations  (27-35).  Everywhere  in  region  I (the  upper  half-space), 
the  reflected  field  is: 


Er  = eyy  l S J fq  RC  sin(fx)  e"gz  df 
x m k^  + q 


(38) 


e;  = - eYy  I (r*j2m./s)Sm  J 3 »C  cos(f»|  e'*  df 


(39) 


E^  = eyy  l Sm  / RC  cos(fx)e"gz  df 


(40) 


-j  nl00  (cosh  gh  - cos  khl 
^ it  gs  sin  kh 


(41) 


Everywhere  in  region  II  fthe  lower  half-space),  the  transmitted  field 
Is 


E1  = . eyy  l S / f9  TC  sin(fx)eqz  df 

q - y2 


(42) 
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Figure  2.  Periodic  array  of  vertical  dipoles  over  the  flat  earth. 


E*  = eYy  l <y*J2w/s)  Sm  S 3_[C  cosffxle^  if 
y q ' y2 

(43) 

Ez  “ ^ I sm  / TC  cos(fx)eqz  df 

(44) 

q2  a - (y  + J 2nm/s)2  + f2 

(45) 

2 2 

Y2  * jo)M2o2  - w m2e2  . 

(46) 
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It  is  found  that  h£  = 0 and  = 0.  The  tangential  components  of  E 
and  H must  be  continuous  across  the  air-earth  Interface.  By  enforcing 
these  boundary  conditions,  we  determine  the  reflection  coefficient  R 
and  the  transmission  coefficient  T: 


(gec  - qc0)e-2gd 
gec  + geo 


(47) 


T = 


2 ge  e gc*  eg^ 
o 

gec  + geo 


(48) 


ec  = eg  ~3c2^  • 


(49) 


A 


Now  that  expressions  are  available  for  the  field  reflected  from 
the  air-earth  interface,  we  are  in  a position  to  evaluate  the  term  z£ 
in  the  transcendental  equation. 

The  reflected  field  makes  the  following  contribution  to  the  voltage 
V(  ) in  Equations  (7)  and  (13): 

Zo — / sin  k ( h- 1 z | )e!"(0,0,z)  dz  . (50) 

0 I sin  kh  -h  z 

oo 

From  Equations  (40)  and  (50), 


Zo  = T0"2sin  kh 


(51) 


Specific  numerical  results  derived  from  the  above  theory  are  given 
in  Section  I IB. 
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2.  The  Finite  Line 


The  foregoing  theory  is  suited  for  ideally  periodic  structures 

i.e.,  those  of  infinite  extent.  If,  however,  we  wish  to  investigate 

the  effects  of  bends  in  the  line,  or  of  various  terminations  and  feeds, 

or  of  an  obstacle  placed  near  the  line,  it  is  appropriate  to  apply  the 
41-42 

method  of  moments  to  a Yagi  line  of  finite  extent.  By  the  method 
of  moments  we  mean  nothing  more  than  a computer-oriented  technique  where 
by  the  integral  equation,  derived  from  Maxwell's  equations  subject  to 
the  appropriate  boundary  conditions  at  the  dipoles,  is  approximated 
by  a set  of  simultaneous  linear  algebraic  equations  whose  unknowns  are 
the  currents  induced  on  the  dipoles  by  some  given  primary  (voltage) 
source.  In  this  method,  each  dipole  is  considered  to  be  divided  into 
N segments  and  a piecewise-sinusoidal  current  of  unknown  amplitude  is 
assumed  on  each  segment.  These  current  amplitudes  form  the  unknowns. 
Written  in  matrix  form,  the  problem  reduces  to  one  of  inverting  the 
kernal  matrix  f usually  an  impedance  matrix)  by  some  algorithm  suited 
to  the  digital  machine,  whereby  one  obtains  numerical  values  for  the 
current  distributions  on  each  dipole. 

Under  this  contract  we  considered  by  the  above  technique  only  Yagi 
lines  in  free-space;  the  earth  can  be  incorporated  if  one  is  willing 
to  accept  considerable  complications  in  the  formulation  of  the  fields. 

Up  to  200  dipoles  were  considered,  but  most  computations  were  made  with 
fewer  than  100  dipoles  and  in  some  cases  as  few  as  30  dipoles. 

The  primary  source  usually  consisted  of  a voltage  source  at  the 
center  of  a dipole,  or  several  sources  in  adjacent  dipoles  phased  to 
correspond  to  the  phase  velocity  of  the  propagating  surface  wave. 

Specific  numerical  results  involving  bends  in  the  Yagi  line,  in- 
truder effects  and  feed  and  termination  designs  are  given  in  Section 
I IB . 
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B.  Some  Calculated  Results 


Because  the  Yagi  line  has  so  many  parameters  available  to  the  de- 
signer (one  of  the  reasons  we  selected  this  structure)  it  was  difficult 
under  this  contract  alone  to  amass  all  the  data  necessary  to  form  a 
complete  picture  of  the  line's  behavior.  In  this  section  we  present 
certain  calculated  data,  based  upon  the  approaches  just  described,  to 
give  some  idea  of  the  line's  potential  as  an  intruder  sensor.  Some 
of  the  computer  codes  employed  are  given  in  References  1 and  2,  should 
more  data  be  desired;  other  computer  codes  (specifically,  those  involv- 
ing the  finite  line)  are  not  documented  explicitly. 

Figures  3 and  4 are  contour  plots  of  the  relative  phase  velocity 
and  attenuation  constant  of  the  primary  mode  supported  by  the  infinitely 
long  Yagi  line  of  copper  rods  of  radius  0.001  operating  at  300  MHz 
in  free  space.  Figures  5-8  present  the  same  data  in  a different  form. 
Comparing  Figures  5-8  with  Figures  9-12  will  give  some  feeling  for  the 
effects  of  reducing  the  dipole  radii  from  0.001X  to  0.00067X,  all  other 
parameters  being  fixed. 

All  the  above  figures  show  the  same  trends:  (1)  the  surface  wave 
slows  down  with  increasing  dipole  length  and  with  decreasing  dipole 
spacing,  all  other  parameters  being  fixed,  (2)  slowing  down  the  wave 
either  by  increasing  dipole  length  or  by  decreasing  dipole  spacing  in- 
creases the  attenuation  constant  due  to  increased  copper  loss,  (3)  slow- 
ing down  the  wave  binds  the  surface  wave  more  tightly  to  the  structure. 
The  last  fact  may  be  deduced  from  Equations  (27 ) - (37 ) and  based  upon 
these  equations  we  have  calculated  and  plotted  fields  for  a few  cases. 
Figures  13-18  present  magnitude  of  the  electric  field  component  Ez  par- 
allel to  the  dipoles  as  a function  of  the  transverse  distance  from  the 
dipoles  along  the  line  y=0,  z=0  (l.e. , in  the  plane  of  a dipole).  The 
dipoles  are  copper  rods  of  radius  0.001X,  the  frequency  is  300  MHz, 
and  the  four  curves  on  each  figure  correspond  to  dipole  lengths  of  2h/x 
= 0.36,  0.40,  0.42  and  0.44.  Dipole  spacings  considered  are  s/X  * 0.30, 
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Figure  4.  Contours  of  constant  attenuation  for  surface  wave  on 
periodic  array  of  dipoles  in  free  space. 
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Figure  5.  Attenuation  constant  versus  dipole  length  for  copper  dipoles 
with  a = 0.001  wavelength,  f * 300  MHz. 
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Figure  8.  Phase  velocity  versus  dipole  length  for  copper  dipoles 
with  a = 0.001  wavelength,  f = 300  MHz. 
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Figure  10.  Phase  velocity  versus  dipole  length  for  copper  dipoles 
with  a = 0.000667  wavelength,  f = 300  MHz. 
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Figure  12.  Phase  velocity  versus  dipole  length  for  copper  diodes 
with  a = 0.00067  wavelength,  f = 300  MHz. 
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Figure  13.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a A =0.001,  sA=0.30,  f=300  MHz. 
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Figure  14.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a/X=0.001,  s/X=0.26,  f=300  MHz. 
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Figure  15.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a A =0.001,  sA=0.-22,  f=300  MHz. 
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Figure  16.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a/A=0.001,  s/A=0.18,  f=300  MHz. 
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Figure  17.  z-coniponent  of  E-field  in  a transverse  cut  in  plane  of  a dipole 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
aA =0.001,  s A =0.14,  f=300  MHz. 
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Figure  18.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp 
a/*=n.om,  sA=ojn,  f=3oo  MHz. 
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0. 26,  0.22,  0.18,  0.14,  and  0.10.  Each  field  plot  is  normalized  to 
correspond  to  1 ampere  of  current  flowing  at  the  center  of  the  reference 
dipole.  Although  this  is  a convenient  normalization  it  is  not  a par- 
ticularly good  one  for  design  purposes.  Better  would  be  a normalization 
such  that  each  field  plot  corresponds  to  a surface  wave  conveying  one 
watt  of  power  across  the  transverse  plane  y = 0 in  which  the  fields  are 
calculated.  By  properly  combining  Equations  (32)  and  (35)  of  Reference 

1 , one  may  derive  the  relationship, 


where  Pt  is  the  power  transferred  across  the  transverse  plane  y=0 

P,  - |l/  R,  iM-sHWcoslkh)]  is  the  power  d1ss1psted  1n 

2tt  ka  sin^(kh)  the  dipole  located  at  y=0  and 

carrying  a current  I at  its 
center  and  exhibiting  a surface 
impedance  Z whose  real  part 
is  R$.  5 

If  R$  is  calculated  for  each  dipole  of  interest  (it  involves  Kelvin 
functions  of  the  dipole  radius,  a,  the  frequency,  and  the  conductivity 
of  copper,  o=57  x 10^  Siemens/meter),  and  if  IQ=1  ampere.  Equation  (52) 
leads  to  Table  I for  P^.  in  wa*  s.  Thus,  if  the  reader  wishes,  he  can 
renormalize  Figures  13-18  to  correspond  to  one  watt  transmitted  simply 
by  dividing  their  vertical  scales  by  the  appropriate  numbers  given  in 
Table  I.  He  can  also  obtain  the  attenuation  and  relative  phase  velocity 
for  each  case  depicted. 

Figures  19-22  present  field  plots  similar  to  the  ones  discussed 
above,  normalized  to  1 ampere  flowing  at  the  center  of  the  reference 
dipole,  but  at  a frequency  of  100  MHz.  The  dipoles  were  of  copper  (field 
calculations  for  aluminum  dipoles  are  not  Included  because  they  differ 
very  slightly  from  that  for  copper)  and  of  two  diameters,  2.05  cm  and 
4.08  cm  (corresponding  to  s/X  = 0.0034  and  0.0068,  respectively).  The 
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Table  I 

Power  transferred,  power  dissipated,  attenuation  and  relative 
velocity  on  a free  space  copper  Yagi  line  at  300  MHz.  Reference 
current  of  1 Ampere,  a/X=0.001 


l/\ 

s/X 

Pt(watts) 

v/c 

560 

. 100 

4 060. 

.11  .000130 

. 9862 

560 

. 140 

623H. 

.11  . 000059 
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.11  . 000025 

.9984 

160 

.250 

23 172. 

.11  . 0000 1 0 
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060 
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EZ  MAGNITUDE  (V/M) 

50.  100.  150.  200. 


SM=  0.900  T = 0.0000 

AM=  0.01023  Z»  0.00 

FREQ  (MHZ) = 100.00 

Figure  19.  z-component  of  E-fleld  In  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a/X*0.00341,  s/X-0.3,  f-100  MHz. 
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lv/« 


'■so-  — 


SM=  0.600  r=  0.0000 

flM=  0.01023  1-  0.00 

FREQ  (MHZ) = 100.00 

Figure  ?0.  z-component  of  E-f lei d in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  Is  1 Amp. 
a/X=0.00341,  s/X=0.?,  f=100  MHz. 
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EZ  MAGNITUDE  (V/M) 

100.  150.  200. 


SM=  0.900  Y«  0.0000 

AM=  0.020H6  Z«  0.Q0 

FREQ  (MHZ) * 100.00 


Figure  21.  z-component  of  E-fleld  In  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  Is  1 Amp. 
a/X =0.00682,  $/X*0.3,  f»100  MHz. 
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EZ  MRGNITUDE  (V/M) 

100.  150.  200. 


SM=  0.600  Y=  0.0000 

RM=  0.02046  Z=  0.00 

FREQ  (MHZ) = 100.00 

Figure  ?2.  z-component  of  E-field  in  a transverse  cut  in  plane  of  a dipole, 
normalized  such  that  current  at  center  of  each  dipole  is  1 Amp. 
a/X=0.00fi82,  s/X=0.2,  f=100  MHz. 
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spacings  between  dipoles  were  0.6  m and  0.9  m (corresponding  to  s/x  = 

0.2  and  0.3,  respectively).  The  corresponding  values  for  Pt,  P^,  a, 
and  v/c  are  given  in  Table  II  which,  as  before,  may  be  used  for  renor- 
malization purposes. 

The  Yagi  line  design  parameters  (spacings,  diameters,  lengths) 
selected  for  Figures  21  and  22  correspond  at  100  MHz  to  a 25.98:1  scaled 
version  constructed  of  copper  for  experimental  testing  at  ~2.5  GHz. 

See  Section  IIC  of  this  report  for  details.  We  calculated  fields  cor- 
responding to  Figures  21  and  22  for  the  scaled  version  at  f = 2.598 
GHz  and  found  them  to  be  essentially  unchanged  (allowing  for  the  scale 
factor)  from  those  at  f * 100  MHz. 

All  the  field  plots  of  Ez  presented  in  this  report  were  made  in 

a transverse  plane  containing  a dipole  (y=0).  Similar  plots  for  planes 
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between  dipoles  display  very  little  differences  from  those  shown  . 

(This  is  not  true  of  other  field  components,  however.) 

We  now  turn  our  attention  to  the  Yagi  transmisison  line  above  a 
flat  air-earth  Interface. 

Figure  23  Is  a contour  plot  of  the  attenuation  constant  of  the 
primary  mode  supported  by  the  Infinitely  long  Yagi  line  of  vertical 
copper  dipoles  of  radius  0.001  suspended  one  wavelength  (measured  from 
centerline)  over  the  flat  earth.  The  earth  is  specified  by  the  con- 
stitutive parameters*  e?s2eQ,  u2=w0»  cj=0.01  mho/m.  This  plot  should 
be  compared  with  Figure  4 for  the  same  Yagi  line  in  free  space.  For 
the  line  in  free  space  it  was  noted  that  the  attenuation  constant  de- 
creases as  dipole  length  decreases.  As  the  dipole  length  decreases, 
the  phase  velocity  Increases  and  the  surface-wave  becomes  more  loosely 
bound  to  the  line,  weaker  currents  flow  on  the  dipoles  and  less  power 
Is  dissipated  in  each  dipole  when  the  power  transmitted  via  the  surface- 


*Some  data  on  Earth  media  electrical  characteristics  are  given  In  Ref- 
erence 44. 
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Table  II 

Power  transferred,  power  dissipated,  attenuation  and  relative 
velocity  on  free  space  copper  Yagi  lines  at  100  MHz.  Reference 
current  of  1 Ampere,  aA  =0.0034  and  0.0068. 

a/X=.341000E  -2 


l/\ 

s/X 

P^f watts) 

P^fwatts) 

v/c 

. 3 60 

.200 

2700. 

.02 

. 000005 

.9945 

.3  60 

. 3 00 

ft  00  1 . 

.02 

. 00000! 

.9992 

.4  00 

.2  00 

3157. 

.02 

.000049 

.9509 

.400 

.300 

494. 

.02 

. 000021 

. 9ft  3ft 

.420 

.200 

1 5 1 . 

.02 

. OOOI24 

.0999 

.420 

. 3 00 

14ft. 

.02 

.000084 

. 949ft 

.4  40 

.2  00 

60. 

.02 

. 000335 

.0029 

.440 

.300 

51. 

.02 

. 0002  66 

.0592 

a/X=.682000E  -2 

.3  60 

.200 

949. 

.01 

. OOOOPft 

. 9853 

.360 

.300 

1 (<40. 

.01 

. 000003 

.9964 

.4  00 

.200 

172. 

.0! 

. 000050 

.92  14 

.4  00 

.300 

192. 

.01 

. 000030 

.9633 

.420 

.2  00 

92. 

.01 

.0001  02 

.8457 

.420 

.300 

69. 

.01 

.000091 

. 9042 

.440 

.200 

44. 

.01 

. 000731 

. 72  7ft 

.440 

.300 

Ift. 

.01 

. 000377 

.78  12 

41. 


FREQUENCY : 300  MHz 
DIPOLE  RADIUS:  a » 0.001  X 
DIPOLE  MATERIAL:  SOLID  COPPER 
HEIGHT  ABOVE_TJjE  EARTH:  d - X 
a • 0.010 


a • 0.020 


0.441- 


0.004x> 


0.006^^' 


0.002 -A/ 


yf  a -0.017  NEPERS /m 

PARAMETERS  OF  THE  EARTH: 
PERMITTIVITY:  «2  ■ 2 «0 
PERMEABILITY:  /i2  “ Mo 
CONDUCTIVITY:  <r,  - 0.01  mho/m 


DIPOLE  SPACING  */X 


Figure  23.  Contours  of  constant  attenuation  for  surface  wave 
on  periodic  array  of  vertical  dipoles  over  the 
flat  earth. 


wave  is  held  fixed.  For  dipole  lengths  receding  from  A/2,  this  same 
type  of  behavior  is  observed  in  Figure  23  for  the  Yagi  line  over  the 
flat  earth.  But  here,  as  the  dipole  length  decreases  further,  the  at- 
tenuation constant  increases  again.  The  surface-wave  becomes  more  loosely 
bound,  the  surface-wave  field  penetrates  more  strongly  into  the  earth 
and  more  power  is  dissipated  there.  Thus,  for  each  value  of  spacing 
s and  height  d,  there  is  an  optimum  dipole  length  that  minimizes  a. 

As  one  increases  the  height  d above  the  earth,  the  attenuation  con- 
stant decreases  rapidly  and  the  surface-wave  characteristics  approach 
those  of  the  array  in  free  space. 

We  have  not  constructed  a contour  plot  of  relative  phase  velocity 
for  the  Yagi  line  over  flat  earth.  Evidently  It  would  be  almost  identical 
to  that  shown  in  Figure  3 for  the  free  space  case,  differing  significantly 
only  in  the  region  where  the  phase  velocity  is  nearly  equal  to  the  velocity 
of  light. 

Figures  24-26  show  the  Influence  on  a and  v/c  of  raising  the  Yagi 
line.  Figures  27-29  present  similar  data  which  show  the  effect  of  de- 
creasing the  dipole  radius  from  0.001A  to  0.0067A. 

Figures  30-33  present  field  plots  of  Ez  of  the  same  copper  Yagi 
pertinent  to  Figures  21  and  22  but  spaced  varying  distances  above  a 
lossy  earth  interface.  The  frequency  is  100  MHz,  normalization  in  all 
cases  is  to  one  ampere  at  the  center  of  the  local  dipole,  and  all  plots 
are  in  the  plane  of  the  dipole  (y=0).  Rather  than  plotting  Ez  as  a 
function  of  transverse  distance  x along  a line  (y=0,  z=0)  as  was  done 
in  Figures  21  and  22  for  the  free  space  case,  we  here  have  plotted  E? 
as  a function  of  z along  three  lines,  (y=0,  x=0.375  m),  (y=0,  x=0.75  m) 
and  (y=0,  x=1.5  m)  extending  from  the  ground  upward.  Figure  30  Is  for 
the  center  line-to-ground  distance  d*2.4  m,  which  is  essentially  unchanged 
from  the  free  space  case.  Figures  31-33  are  for  progressively  smaller 


ATTENUATION  (OB/KM) 

o.o  20.0  no.o  ao.o  «o.o  voo.o 


i 1 1 i i 

o.  >a  o.»«  0.40  « 0.44  °-4* 

LENGTH  OF  DIPOLE 

Figure  24.  Attenuation  constant  versus  dipole  length  for  copper  dipoles 
over  earth  (u2=M0,e2«2eo,o2*0.01  mho/m).  f=300  MHz. 
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Figure  26.  Phase  velocity  versus  dipole  length  for  copper  dipoles 
over  earth  (p1=po,e2=2eo,o2=0.01  mho/m).  f=300  MHz. 
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Figure  27.  Attenuation  constant  versus  dipole  length  for  copper  dipoles 
over  earth  fp?=po,e?*2eo,a2=0.01  mho/m).  f=300  MHz. 
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Figure  29.  Phase  velocity  versus  dipole  length  for  copper  dipoles 
over  earth  fy1-uo,e?=2e(),a?=0.01  mho/m).  f=300  MHz. 
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Figure  30.  z-component  of  E-fleld  for  a copper  Yagi  over  lossy  ground 
with  £,=16,  o2=  0.02,  aA =0.00682,  sA=0.3,  «-A=0.42, 
f=100  MHz,  d=2.4  m.  Plots  are  along  lines  parallel  to 
a dipole  at  distances  x=0.375,  0.750,  and  1.5  meters. 
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Figure  31.  z-component  of  E-field  for  a copper  Yagi  over  lossy  ground 
with  £,=16,  o2=0.02,  a/X=0. 00682,  s/X=0.3,  H/X=0.42, 
f=100  MHz,  d=2.0  m.  Plots  are  along  lines  parallel  to 
a dipole  at  distances  x=0.375,  0.750,  and  1.5  meters. 
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Figure  32.  z-component  of  E-field  for  a copper  Yagi  over  lossy  ground 
with  ^=16,  ct2=0.02,  a/\=0. 00682,  s/x=0.3,  Jl/X=0.42, 
f=100  MHz,  d=i.5  m.  Plots  are  along  lines  parallel  to 
a dipole  at  distances  x=0.375,  0.750,  and  1.5  meters. 
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SM-  0.90  Y-  0.000 
RM-  0.02046  X-  0.375 
DM-  1.00  0.750 
FREQ  (MHZ)-  100.0  1.500 


Figure  33.  z-component  of  E-field  for  a copper  Yagi  over  lossy  ground 
with  eo=16,  o2=0.0?,  a/X=0.00682,  s/x=0.3,  S,/X=0.4?, 
f=100  MHz,  d=f.O  m.  Plots  are  along  lines  parallel  to 
a dipole  at  distances  x=0.375,  0.750,  and  1.5  meters. 
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values  of  d.  The  smallest  value  of  d shown  is  one  where  the  lower  tip 
of  a dipole  is  about  X/10  above  the  interface;  smaller  values  of  d are 
felt  to  lead  to  incorrect  results  due  to  the  fact  that  our  model  does 
not  allow  skewing  of  the  current  distribution  on  the  dipole  by  the  pres- 
ence of  the  earth.  Values  of  v/c  and  ot  are  tabulated  in  Table  III. 


Table  III 

Calculated  values  of  v/c  and  a for  the  Yagi  lines 
over  ground  described  in  Figures  30-33. 


d(m) 

a ^ km  ) 

c 

00 

0.79026 

0.90417 

2.40 

28.950 

0.90292 

2.00 

65.104 

0.90107 

1.50 

175.34 

0.89402 

1.00 

453.58 

0.86771 

The  presence  of  the  earth  distorts  the  Ez  component,  increasing 
its  value  near  the  earth  relative  to  that  along  the  centerline  (y=0, 
z=0).  This  is  good  from  the  standpoint  of  providing  detection  of  an 
intruder  trying  to  roll  beneath  the  line,  but  bad  from  the  standpoint 
of  attenuation  due  to  earth  losses.  Unfortunately  we  have  no  calculations 
of  earth  losses  comparable  to  of  Table  II,  and  therefore  have  not 
provided  numbers  to  renormalize  Figures  30-38  to  one  watt  of  transferred 
power. 

Another  aspect  of  the  Yagi  line  (or  of  any  open  transmission  line 
for  that  matter)  is  the  effect  of  a curve  upon  the  surface  wave  propa- 
gating along  the  line.  Invariably  the  discontinuity  caused  by  the  curve 
will  cause  a reflection  which,  together  with  radiation  of  energy  as 
the  wave  proceeds  around  the  curve,  will  result  in  a reduction  of  energy 
transmitted  by  the  surface  wave  propagating  beyond  the  curve.  A high 
transmission  loss  is  undesirable  and  it  is  of  interest  to  know  if  and 
how  one  can  design  a line  to  convey  the  surface  wave  around  a curve 
with  tolerable  loss.  We  had  a tentative  look  at  this  problem. 
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Figure  34  represents  a top  view  of  a Yagi  line  of  90  dipoles  in 
free  space,  dipoles  number  1-30  and  61-90  forming  straight  input  and 
output  sections,  respectively,  the  remainder  forming  a 90°  circular  arc. 

The  dipoles  all  have  a spacing  S=0.3X,  length  d=0.44X  and  a radius  a=0.001x. 
Dipole  No.  1 is  excited  with  1 ampere,  causing,  as  described  previously, 
surface  waves  and  "strange"  waves  to  travel  from  the  source  to  the  end 
of  the  line  (outgoing  waves)  and  from  the  end  of  the  line  back  to  the 
source  (reflected  waves).  By  solving  for  all  the  currents  on  the  90 
dipoles  by  the  method  of  moments  technique  we  can  choose  sets  of  6 dipoles, 
first  in  the  straight  input  section,  then  in  the  curve,  and  finally 
in  the  straight  output  section.  Each  of  these  will  yield  the  amplitudes 
of  the  currents  and  propagation  constants  of  the  surface  waves  and  strange 
waves  proceeding  in  both  directions  in  each  of  these  three  regions. 

Figure  35  shows  the  results  for  the  outgoing  and  reflected  surface  waves 
only  (the  outgoing  and  reflected  strange  waves  are,  respectively,  one 
and  two  orders  of  magnitude  below  the  reflected  surface  wave  and  are 
not  presented  here).  The  solid  curve  CQ  represents  the  (eye- Interpolated) 
curve  passing  through  the  current  amplitudes  associated  with  the  out- 
going surface  wave.  Clearly,  these  currents  experience  a substantial 
drop  in  their  procession  around  the  curve,  indicating  substantial  trans- 
mission loss  due  to  radiation  (wave  attenuation  due  to  copper  loss  is 
negligible  by  comparison).  The  dashed  curve  Cr  represents  the  (eye- 
interpolated)  curve  passing  through  the  current  amplitudes  associated 
with  the  reflected  surface  wave;  an  identical  percentage  transmission 
loss  is  experienced  by  this  wave  as  it  proceeds  around  the  curve,  as 
would  be  expected  from  the  symmetric  geometry  of  the  whole  line.  The 
two  curves  marked  C + C and  C - C are  simply  the  sum  and  difference 
of  CQ  and  Cr.  The  heavy  dots  represent  the  actual  current  amplitudes 
as  obtained  by  the  method-of -moments  before  resolution  Into  surface 
and  strange  modes.  The  fact  that  a few  points  fall  outside  the  region 
bounded  by  CQ  + Cr  and  C0  + Cr  is  evidence  that  strange  modes  exist 
but  are  small.  It  may  be  commented  here  that  the  velocities  of  the 
strange  modes  were  found  to  exceed  the  speed  of  light  - as  would  be 
expected  if  they  were  due  to  direct  radiations  from  the  ends  of  the  line. 


55 


Figure  35.  Transmission  around  a 90°  arc  of  a finite  Yagi  line. 


that  radiation  proceeding  in  straight  lines  to  all  points  in  the  curved 
Yagi  at  the  speed  of  light.  Then  testing,  as  we  did,  along  a line  skew 
to  these  directions  of  radiation,  wavelengths  larger  than  actually  present 
would  be  encountered  which  would,  in  our  model,  be  interpreted  as  strange 
waves  proceeding  along  the  Yagi  at  speeds  greater  than  that  of  light. 

In  fact,  this  observation  lends  credence  to  our  speculation  that  the 
strange  waves  are  indeed  due  to  radiation  from  the  feed  region  and  re- 
flection region. 

In  Figure  35  all  dipole  spacings  are  held  fixed.  We  speculated 
that  if  the  surface  wave  guided  by  the  Yagi  line  were  more  tightly  bound 
to  the  structure  in  the  vicinity  of  the  curve,  less  would  be  lost  to 
radiation  and  transmission  loss  would  be  reduced.  To  test  this  specu- 
lation, the  spacings  of  the  dipoles  in  the  90°  curved  region  were  re- 
duced from  S=0.3X  to  S=0.15X  in  increments  of  0.05X  (the  smaller  spac- 
ings require  more  of  the  dipoles  to  lie  in  the  curved  region,  leaving 
fewer  of  the  total  number  of  dipoles  to  lie  in  the  straight  regions). 
Figure  36  shows  the  interpolated  dipole  current  magnitudes  (in  the  curved 
region  only)  for  the  outgoing  surface  wave.  There  is  clearly  reduced 
transmission  loss  as  the  dipole  spacings  are  reduced  (a  similar  effect 
would  be  expected  if  the  spacings  were  held  fixed  and  the  dipoles  in- 
creased in  length,  the  essential  feature  being  a more  tightly  bound 
wave  in  the  vicinity  of  the  curve).  Thus,  we  feel  confident  that  the 
Yagi  line  inherently  has  design  parameters  which  allow  substantial  con- 
trol of  the  field  structure  to  accommodate  curves,  etc.  Whether  or 
not  these  are  compatible  with  the  aim  of  being  a good  intruder  sensor, 
we  have  not  gathered  enough  information  to  know. 

The  last  two  figures.  Figures  37  and  38,  show  the  (negligible) 
effect  of  a linearly  tapered  spacing  of  dipoles  in  the  input  and  out- 
put regions  adjacent  to  the  curve.  The  thought  here  was  that  the  dis- 
continuity caused  by  the  curve  and  sudden  reduction  in  dipole  spacing 
might  cause  a non-negllgible  reflection  of  the  outgoing  surface  wave 
as  it  enters  the  curved  region.  The  tapers  were  intended  to  mitigate 
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Figure  37.  Transmission  around  a 90  arc  of  a finite  Yagi  line 


this  reflection.  Figure  38  (currents  plotted  vs  the  angle  <J>  around 
the  curve  rather  than  vs  dipole  number  as  in  Figure  36)  indicates  very 
little  improvement  due  to  the  taper. 

The  method  of  moments  technique  can  be  applied  to  a finite  Yagi 
line  in  the  presence  of  an  "intruder"  to  estimate  the  reflection  that 
would  be  caused  by  the  intruder.  Figures  39  and  40  illustrate  typical 
situations.  A finite  Yagi  line  in  free  space  is  postulated  which  is 
excited  in  a chosen  manner  - in  Figures  39  by  a single  center-fed  dipole, 
or  in  Figure  40  by  10  dipoles  center  fed  by  voltage  sources  chosen  such 
that  the  currents  they  induce  on  their  respective  dipoles  are  constant 
from  one  excited  dipole  to  the  next  and  exhibit  a progressive  phase 
corresponding  to  the  velocity  of  the  surface  wave  supported  at  the  ex- 
citation frequency.  To  reduce  any  reflections  at  the  source  end,  8 di- 
poles are  center-terminated  in  resistances,  increasing  from  10  ohms 
at  the  dipole  nearest  the  excited  region  to  80  ohms  at  the  farthest 
dipole.  To  reduce  reflection  at  the  far  end,  a similarly  tapered  lossy 
load  arrangement  is  used.  In  the  absence  of  an  intruder,  th~  current 
IQe~^y  on  each  dipole  in  the  mid  section  of  the  line  corresponds  to 
the  appropriate  surface  wave  proceeding  from  source  to  load.  As  in- 
dicated by  the  dashed  line  in  Figure  39,  this  incident  wave  is  quite 
free  of  a reflected  component,  verifying  that  our  tapered  load  was  doing 
a very  good  job  of  terminating  the  surface  wave.  An  intruder  (in  the 
form  of  a thin  perfectly  conducting  wire  introduced  at  dipole  #20  in 
Figure  39  and  dipole  #40  in  Figure  40)  causes  reflected  surface  wave 
currents  to  be  induced  on  the  Yagi,  evidenced  by  the  standing  wave  pat- 
tern shown  as  a solid  line  in  Figure  39.  The  dots  in  this  curve  are 
the  values  of  current  magnitudes  as  obtained  from  the  method  of  moments 
while  the  solid  line  is  merely  a hand-drawn  interpolation  of  the  standing 
wave.  A more  careful  study  was  made  of  the  reflection  coefficient  in 
the  configuration  of  Figure  40.  Here,  only  the  currents  on  the  midsection 
of  the  lines  (dipoles  25-35)  are  presented  because  it  is  in  this  region 
that  we  expect  to  see  a combination  of  Incident  and  reflected  waves. 
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Figure  40.  Currents  induced  on  a finite  copper  Yagi  line  in  free 
space  in  the  presence  of  a nearby  intruder.  The 
intruder  is  a dipole  of  height  0.44A  at  a distance 
x=0.2A  from  the  40th  dipole  in  the  line.  f=300  MHz, 
line  parameters  given  in  Fig.  39. 
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The  dots  are  the  values  of  the  currents  as  obtained  from  the  method  of 
moments,  the  solid  line  is  the  curve  IQe'Yy  + RIQeYy,  where  the  values 
of  the  two  unknowns,  I and  R,  were  obtained  by  matching  this  mathe- 
matical form  to  the  actual  currents  obtained  by  method  of  moments.  This 
matching  was  performed  in  a least-mean  square  sense  over  the  11  dipoles 
indicated  in  Figure  40.  Comparison  of  the  dots  with  the  solid  curve 
lead  us  to  believe  that  our  representation  of  the  currents  as  incident 
and  reflected  surface  waves  with  reflection  coefficient  magnitude  |R| 
~0.14  is  reasonably  accurate. 

With  the  foregoing  comments  in  mind.  Figures  41-44  present  the 
magnitude  of  the  surface-wave  reflection  coefficient  of  an  intruder  near 
an  aluminum  Yagi  in  free  space  at  100  MHz.  Two  dipole  radii,  a,  and 
two  dipole  spacings,  s,  are  considered  for  a variety  of  dipole  lengths, 
l.  The  intruder  in  all  cases  is  a thin  perfectly  conducting  rod  of 
length  0.44A  (1.3?  meters)  which  is  parallel  with  and  symmetrical ly 
located  with  respect  to  the  Yagi  dipoles  and  up  to  A/2  (1.5  meters) 
away  from  them.  The  reflection  coefficient  in  all  cases  is  referenced 
to  the  Yagi  dipole  nearest  the  intruder.  The  curves  appear  to  flatten 
out  at  distances  larger  than  about  0.3A,  particularly  for  those  com- 
binations of  s and  l (smaller  s,  larger  £)  which  correspond  to  slower 
surface  waves,  i.e.,  more  tightly  bound  waves.  This  indicates  that 
we  are  observing  the  effects  of  direct  radiation  from  the  feed  region 
to  the  intruder  region  which  corrupts  the  assumptions  upon  which  our 
definition  of  reflection  coefficient  is  based.  The  results  are  probably 
reliable  for  values  of  x/A  <_  0.3,  i.e.,  within  1 meter  of  the  line. 

An  additional  comment  is  in  order  regarding  the  computational  in- 
vestigation just  described,  and  that  involves  the  radiation  efficiency 
of  the  entire  configuration  under  consideration.  By  radiation  efficiency 
we  mean  the  ratio  of  power  radiated  to  total  power  supplied  by  the  pri- 
mary feed.  Radiation  takes  place  from  the  feed  region  and  termination 
region  and,  where  an  intruder  is  present,  from  the  intruder  and  intruder 
region.  Figure  39,  for  example  shows  that  59. IX  of  the  input  power 


Figure  42.  Reflection  coefficient  of  the  surface  wave  on  Yagi  lines  due 
to  an  intrude1"  (wire)  as  a function  of  intruder  distance. 


Figure 48s  Reflection  coefficient  of  the  surface  wave  on  Yagi  lines  due 
to  an  intruder  (wire)  as  a function  of  intruder  distance. 


Reflection  coefficient  of  the  surface  wave  on  Yagi  lines  due 
to  an  intruder  (wire)  as  a function  of  intruder  distance. 


went  into  radiation  and  Introduction  of  the  intruder  raised  that  to 
63. 7%.  Much  of  this  radiation  is  due  to  the  feed  configuration,  in 
this  case  a single  dipole  backed  by  a linearly  tapered  resistive  load. 

If  a distributed  feed,  such  as  that  described  previously  for  Figure 
40,  is  used,  we  would  expect  a better  transfer  of  power  into  surface 
wave  fields  and  less  radiation,  i.e.,  a smaller  radiation  efficiency. 
This  criterion,  small  radiation  efficiency,  was,  in  fact,  the  one  used 
to  arrive  at  the  design  consisting  of  a distributed  feed  backed  by  a 
linearly  tapered  load  used  in  Figures  40-44.  In  Figure  45  such  a feed 
together  with  a linearly  tapered  load  is  assumed  at  the  ends  of  a copper 
Yagi  line  in  free  space  whose  length  is  varied  to  include  from  30  to 
40  dipoles.  The  radiation  efficiency  hovers  about  35%.  We  speculate, 
though  we  have  not  investigated  it  computationally,  that  providing  more 
driven  elements  and  including  more  dipoles  in  the  termination  regions 
will  reduce  this  radiation  efficiency  further,  though  never  to  zero 
in  a practical  case. 

C.  Some  Experimental  Results 

To  verify  theoretical  results,  some  experimental  data  were  taken 
on  a scaled  Yagi  transmission  line.  If  the  full  sized  version  is  as- 
sumed to  operate  at  100  MHz,  the  scaled  model  was  designed  to  operate 
at  2.598  GHz,  i.e.,  it  was  25.98  times  smaller.  Choosing  copper  dipoles 
for  both  the  full  sized  line  and  the  scaled  version  disturbs  the  exact 
scaling  insignificantly  because  conductivity  scales  almost  1:1  for  ex- 
cellent conductors. 

The  Yagi  dipoles  were  supported  by  a thin  ribbon  of  polyfoam  (er 
- 1.05)  whose  width  was  about  half  the  dipole  length.  At  the  design 
frequency  of  2.598  GHz  the  dipoles  had  length  0.44X  U=2"),  diameter 
0.0134A  (2a=l/16"),  and  the  dipole  spacing  was  0.3X  (s=1.36").  This 
structure  was  carefully  suspended  on  a system  of  very  fine  1 lb  test 
nylon  string  zig-zagged  across  a trough-like  wooden  structure  2 feet 
wide,  1 foot  deep  and  12  feet  long.  For  convenience,  the  dipoles  were 
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Figure  45.  Radiation  efficiency  of  a finite  copper  Yagi  line  as  a 
function  of  its  length.  f=300  MHz. 
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oriented  horizontally  as  shown  in  Figure  46.  At  one  end  the  line  was 
terminated  in  a 1 ft  x 1.5  ft.  short  circuiting  aluminum  plate  (not 
shown);  at  the  other  end,  the  line  exited  through  an  aperture  in  a wall 
of  8"  thick  hairflex  absorber  (shown  in  the  photograph)  and  under  its 
own  weight  proceeded  to  curve  gently  upward  over  a 5 foot  length  toward 
a log  periodic  antenna  from  which  it  was  suspended,  as  shown  in  Figure 
47.  As  suggested  by  our  calculated  data,  the  dipole  spacing  was  reduced 
to  0.2A  (s=0.91")  in  this  curved  region  to  better  trap  the  surface  wave 
as  it  proceeded  from  the  source  to  the  line  beyond  the  aperture.  The 
log  periodic  antenna  was  simply  one  which  was  available  and  worked  quite 
well.  As  expected  the  direct  radiation  from  the  feed  antenna  was  re- 
duced in  the  measurement  region  by  this  combination  of  geometrical  ori- 
entation and  absorbing  wall  and,  except  for  positions  very  near  the 
aperture,  only  a surface  wave  field  was  detected  on  the  line,  as  de- 
sired. 

A balanced  loop  antenna,  0.5"  in  diameter,  carefully  constructed 
of  very  fine  coaxial  cable,  served  as  a probe  whose  output  was  carried 
vertically  upward  (a  line  of  zero  tangential  E-field  when  centered  over 
the  Yagl ) to  slightly  larger  coaxial  cable,  thence  to  a moveable  car- 
riage and  eventually  to  a Scientific  Atlanta  superheterodyne  receiver. 

As  shown  in  Figure  48,  the  carriage  consisted  of  a pair  of  modified 
model  railroad  flatcars  riding  on  two  sets  of  tracks  along  a wooden 
shelf  which  by  construction  was  very  level.  A micrometer  support  on 
the  carriage  permitted  the  loop  to  be  moved  in  two  dimensions  in  a plane 
transverse  to  the  surface  wave  propagation  direction.  By  carefully 
laying  the  Yagi  line  parallel  to  the  tracks  and  adjusting  the  tautness 
of  the  support  strings  so  that  the  line  maintained  a uniform  distance 
below  the  probe,  a type  of  "slotted  line"  was  formed  for  the  measurement 
of  attenuation  and  wave  velocity. 

When  centered  over  the  Yagi  line  and  oriented  in  a plane  of  a di- 
pole, the  loop  detected  a signal  essentially  proportional  to  the  longi- 
tudinal component  of  H (H^)  and,  being  balanced,  it  was  insensitive 
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Figure  46.  Photograph  of  experimental  Yagi  line. 


Top*- 


Figure  47.  Photograph  of  log  periodic  feed  antenna. 
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Figure  48.  Photograph  of  carriage  and  loop  probe. 


to  the  electric  field.  The  Hy  field  component  varies  greatly  along 
the  line,  being  strongest  aside  the  dipoles  and  quite  weak  between  them. 
Thus,  the  loop  as  oriented  detected  signals  essentially  proportional 
to  the  current  in  a particular  dipole  when  it  was  located  adjacent  to 
that  dipole.  When  moved  either  side  of  that  location,  the  signal  dropped 
precipitously.  In  a sense  then,  ours  was  a sampled  data  system,  where 
the  uniform  spatial  sampling  interval  was  set  by  the  dipole  spacing, 
and  the  voltage  outputs  were  proportional  to  the  amplitudes  of  the  stand- 
ing wave  at  the  respective  dipole  locations  on  the  short  circuited  line. 
Since  the  sampled  standing  wave  has  a period  different  from  the  sampling 
period,  the  recorded  samples  exhibited  the  periodicity  of  a beat  wave 
(much  like  the  FFT  cf  a periodic  signal).  The  period  of  this  beat  wave 
can  be  extracted  by  eye  as  the  envelope  of  the  recorded  data  and  it, 
together  with  knowledge  of  the  dipole  spacing  and  frequency,  permits 
determination  of  the  wavelength  of  the  surface  wave  and  hence  its  phase 
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velocity.  As  an  example  of  raw  data  at  f =2 . 519  GHz,  Figure  49  shows 
a graph  of  bars  whose  amplitudes  are  proportional  to  the  log  of  the 
power  received  at  39  dipole  stations.  The  envelope  of  these  bars  was 
sketched  by  eye  and  found  to  have  an  average  period  of  P=4.5  dipole 
spacings.  It  may  be  shown  that  the  relative  velocity  is  given  by  the 
expression 


v 2s  f 
c " c 1 - I/P  • 


and  for  the  line  we  constructed  this  becomes. 


(53) 


| = 0.2309  x 10'9  T _ [yp  . (54) 

Using  data  such  as  those  shown  in  Figure  49  and  applying  Equation 
(54)  we  arrive  at  the  experimentally  determined  values  of  v/c  displayed 
as  dots  in  Figure  50.  Values  of  v/c  predicted  by  our  mathematical  model 
are  shown  as  the  dashed  curve.  They  show  the  same  trend  but  are  too 
high  by  15-20^.  This  difference  is  at  least  partly  explained  by  the 
fact  that  the  dipoles,  being  rather  thick  (£/2a=32),  act  longer  at  the 
test  frequency  than  the  2"  they  appear  to  be  optically;  being  effectively 
longer,  the  theory  indicates  that,  indeed,  the  surface  wave  should  be 
slower  than  predicted.  The  difficulty  in  the  theory  arises  from  the 
current  distribution.  Equation  19),  assumed  a priori  to  exist  on  each 
dipole.  In  this  simple  model,  the  dipole  is  considered  to  be  segmented 
into  two  equal  parts  of  length  h with  piecewise  sinusoidal  currents 
on  each,  matched  at  the  center.  A more  accurate  description  of  the 
current  would  be  given  by  an  N-segment  model,  in  which  case  the  current 
distribution  is  not  assumed  a priori  but  is  determined  (within  the  reso- 
lution permitted  by  the  N-segment  basis  set)  along  with  the  propagation 
constant.  Because  this  model  better  resolves  the  current  distribution 
near  the  ends  of  the  dipoles,  one  might  expect  it  to  yield  more  accurate 
values  of  v/c  and  a.  The  price  one  pays  for  increased  N is  a more  complicated. 
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Figure  50.  Calculated  and  experimental  values  of  v/c  and  a over  a frequency  band. 


slower  algorithm  for  determining  the  propagation  constant.  As  a com- 
promise alternative,  the  two  segment  model  can  be  retained  and  an  ef- 
fective length,  £+A£  , can  be  used  rather  than  £ for  the  dipole  length. 
This  was  done  for  an  increase  in  length  of  4%  (£=2.08")  and  the  cor- 
responding theoretical  values  of  v/c  plotted  as  the  solid  curve  in  Fig- 
ure 50. 

There  are  some  additional  points  to  be  made  regarding  the  v/c  data 
in  Figure  50.  First  of  all,  the  mathematical  model  using  an  effective 
dipole  length  £=2.08"  (two-  or  six-segment  model)  misbehaved  above  about 
f=2.55  GHz,  that  is,  calculated  data  points  for  v/c  and  a became  ran- 
domly erratic.  The  algorithm  was  having  difficulty  finding  a complex 
root  of  Equation  (13),  one  of  the  reasons  being  the  incorporation  of 
a divergent  series  in  the  formulation  of  Equation  (13).  (See  discussion 
of  this  point  in  section  I I A. ) A second  formulation,  not  discussed 
in  this  report  but  described  in  a forthcoming  paper45,  avoided  this  di- 
vergent series  and  was  found  to  produce  good  values  of  v/c  and  a beyond 
f=2.55  GHz.  If  the  effective  dipole  length  £=2.08"  (4%  increase  over 
physical  length)  is  retained,  this  convergent  series  formulation  with 
N=2  yields  the  dotted  curves  for  v/c  and  a shown  in  Figure  50.  If  the 
effective  length  £=2.086"  is  used  (4.3%  increase  over  physical  length) 
with  N=2,  the  dashed  curves  for  v/c  and  a result.  In  addition,  the 
method-of-moments  technique  was  applied  to  a finite  Yagi  of  the  given 
dimensions  (incorporating  a 4%  increase  in  dipole  length)  as  described 
in  section  I IB  and  values  of  v/c  and  a were  found  consistent  with  those 
given  by  the  dotted  curves  in  Figure  50.  Within  the  approximations 
inherent  in  our  theory  (a  two-segment  piecewise  sinusoidal  dipole  cur- 
rent model,  use  of  an  empirical  effective  length  for  the  dipoles,  no 
effects  of  the  foam  ribbon  accounted  for,  etc.)  we  will  accept  the  dashed 
curve  of  v/c  and  a as  being  satisfactory. 

As  Figure  50  shows,  the  surface  wave  slows  down  as  frequency  in- 
creases. At  that  frequency  f where  the  dipoles  are  a half  guide  wave- 
length long,  i.e.,  where  Bs=tt,  the  currents  on  adjacent  dipoles  are 
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180°  out  of  phase,  which  would  correspond  to  waves  proceeding  in  either 
the  +y  or  -y  direction  along  the  Yagi.  This  nonuniqueness  implies  that 
in  fact  no  wave  propagates  in  the  +y  direction  along  the  line  when  gs 
- it  and  that  the  frequencies  are  cutoff  frequencies  for  the 

surface  wave  mode  in  question.  These  values  of  f are  indicated  by 
the  sloped  line  in  Figure  50;  our  theoretical  data  hugs  this  line  fairly 
well,  with  a little  undershoot. 

The  region  around  the  cutoff  frequency  is  the  first  Bragg  inter- 
action region.  If  the  Yagi  line  were  lossless,  the  surface  wave  would 
abruptly  change  its  character  at  a unique  f from  a propagating  wave 
whose  group  velocity  becomes  zero  at  f to  an  evanescent  wave  which 
does  not  propagate  but  only  attenuates.  Since  loss  must  occur  in  any 
realistic  case,  this  sharp  division  between  "pass"  and  "stop"  bands 
is  blurred  and  this  is  displayed  in  our  theoretical  curves  by  a smooth 
merging  of  the  two  regions. 

Experimental  data  taken  over  the  entire  12  foot  Tine  exhibited 
the  influence  of  loss,  i.e.,  diminishing  standing  wave  ratio  as  one 
proceeded  from  short  circuit  to  source.  In  the  passband  frequency  range, 
f=2.41  to  2.5  GHz,  this  effect  was  too  small  and  the  measurement  in- 
accuracies too  great  to  resolve  experimental  values  of  a with  confidence. 
In  the  stopband  above  f=2.56  GHz,  however,  where  g remains  essentially 
constant  and  attenuation  is  the  dominant  feature  of  the  Yagi  field, 
instead  of  standing  waves,  a severe  attenuation  of  the  field  was  ob- 
served as  one  proceeded  away  from  the  source  end.  Figure  51  shows  the 
recorded  bar  graph  and  approximating  envelope  taken  at  f=2.574  GHz. 
Because  the  vertical  scale  is  in  dB  the  slope  of  the  (approximately) 
linear  envelope  is  a measure  of  the  attenuation  constant  a.  The  crosses 
in  Figure  50  display  a few  experimental  values  obtained  in  this  way 
for  comparison  with  the  aforementioned  theoretical  curves  of  a. 
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Figure  51.  Representation  of  experimental  data  recorded 
for  Yagi  line  at  f=2.574  GHz. 

A limited  amount  of  experimental  data  were  accumulated  for  the 
Yagi  line  placed  near  a lossy  interface,  as  shown  in  Figure  52.  The 
lossy  material  was  hairflex  absorber  with  its  backside  nearer  the  line. 
Typical  recorded  data  in  the  form  of  bar  graphs  are  not  presented:  it 
may  be  said,  however,  that  for  the  high  attenuation  cases,  the  data 
appeared  very  similar  to  Figure  51  with  very  little  evidence  of  a standing 
wave,  making  evaluation  of  v/c  impossible;  for  the  low  attenuation  cases, 
the  data  appeared  more  like  Figure  49,  with  much  standing  wave  evident 
and  attenuation  so  small  that  it  became  difficult  or  impossible  to  ex- 
tract values  of  a.  In  between,  the  recorded  data  showed  an  attenuated 
standing  wave  near  the  short  circuit,  from  which  v/c  could  be  deduced, 
and  a linear  rise  near  the  source  end  from  which  a could  be  extracted. 
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Figure  52.  Photograph  of  experimental  Yagi  line  with  hairflex  inter- 
face in  place. 

Figures  53  and  54  show  experimental  values  of  v/c  and/or  a,  re- 
spectively, obtained  in  this  way  at  three  frequencies  in  the  pass  band 
region  for  three  nominal  centerl ine-to- interface  distances,  d - 3.175 
cm,  3.81  cm  and  10  cm.  The  solid  curve  in  Figure  53  is  the  same  as 
that  in  Figure  50  for  v/c  of  the  free  space  case.  As  expected,  we  notice 
that  the  influence  of  the  lossy  interface  is  primarily  on  a and  very 
little  on  v/c,  especially  for  the  largest  height  d=10  cm.  The  fuzzy 
nature  of  the  interface  made  precise  evaluation  of  distance  d impossible 
and,  in  fact,  this  distance  varied  a bit  from  point-to-point  along  the 
line  about  the  nominal  values  quoted  above.  Also,  we  do  not  know  the 
values  of  dielectric  constant  and  conductivity  of  the  near  side  of  the 
hairflex  panels.  These  two  intangibles  made  any  theoretical  curves 
generated  for  the  Yagi  line  over  the  earth  of  doubtful  value  for  com- 
parison with  experimental  data.  We  did  try,  however.  Using  estimated 
values  for  hairflex  dielectric  constant  and  conductivity  of  e2=2e0, 
xj?=0.01  Siemens/m,  respectively,  we  calculated  v/c  and  a for  a range 
of  heights  about  the  value  d=l%"  (3.81  cm).  Table  IV  displays  some 
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RELATIVE  PHASE  VELOCITY  ( v/c ) 


of  the  results.  Increasing  o^  to  0.02  Siemens/m  increased  the  values 
of  d only  slightly. 


Table  IV 

Comparison  of  theoretical  and  experimental  values  of  v/c 
and  a for  Yagi  over  lossy  interface. 


f (GHz) 

d= 

3.81  cm 

d*5 

cm 

Experimental 
d=3.81  cm 

v/c 

a(— ) 
Vm  ) 

v/c 

v/c 

El 

2.432 

EE9 

mm 

7.3 

0.83 

1 

2.499 

ESS 

1 

EH 

5.4 

0.78 

4.86 

2.550 

0.64 

18.8 

0.66 

3.3 

0.67 

7.1 

III.  THE  LEAKY  COAXIAL  CABLE 

A.  Synopsis  of  Theoretical  Approach 

In  Reference  46  theories  are  presented  for  calculating  the  prop- 
agation constants  and  electromagnetic  field  distributions  of  waves  sup- 
ported by  three  types  of  transmission  lines.  Respectively,  they  are: 
the  conducting  wire  surrounded  by  a dielectric  sleeve47  (G-line),  the 

coaxial  cable  with  center  conductor  periodically  interrupted  by  complete 

40  40 

circumferential  slots  ’ (C-llne),  and  the  coaxial  cable  with  outer 
conductor  periodically  interrupted  by  incomplete  circumferential  slots 
(R-line).  In  all  cases,  the  lines  were  infinitely  long  and  straight 
and  reside  either  in  free  space  or  parallel  to  a free  space-planar  earth 
interface.  Because  both  the  G-  and  C-lines  are  4>-indpendent  geometries, 
their  assumed  field  structures  are  independent,  whereas  that  of  the 
R-line  must  incorporate  higher  order  <j>-dependent  trigonometric  modes. 
Thus,  although  the  R-line  was  of  primary  Interest,  the  G-  and  C-lines 
served  as  stepping  stones  to  the  more  Intricate  R-line  analysis.  Ad- 
ditionally, independently  derived  theory  and  calculations  existed  for 
the  G-line  over  earth^-^  which  could  serve  as  a check  case. 


Since  all  three  transmission  lines  are  essentially  cylindrical 
in  their  geometry,  the  functions  best  suited  to  the  expansion  of  their 
fields  are  the  cylindrical  wave  functions,  which  serve  quite  well  in 
the  free  space  cases.  The  presence  of  a flat  earth  interface,  however, 
complicates  matters  because  the  cylindrical  wave  functions  are  inappro- 
priate forms  for  matching  boundary  conditions  on  planar  surfaces; 
thus,  more  natural  functions  would  be  plane  wave  functions.  As  was 
done  for  the  Yagi  line  the  presence  of  the  flat  earth  is  accommodated 
mathematically,  first  by  solving  for  the  fields  of  the  free  space  con- 
figuration in  terms  of  cylindrical  wave  functions,  then,  by  transforming 
each  of  these  into  continuous  spectra  of  incident  plane  waves  (Fourier 
integral  expansion)  each  of  whose  interaction  with  the  earth  is  described 
compactly  by  familiar  Fresnel  reflection  and  transmission  coefficients 
and,  lastly,  by  forcing  the  sum  of  the  incident  and  reflected  fields 
to  satisfy  a boundary  condition  back  on  the  transmission  line  proper. 
Since  the  transformations  of  cylindrical  wave  functions  into  plane  wave 
functions  are  common  to  all  three  transmisison  line  analysis  and  are 
not  all  found  in  the  literature,  their  derivations  are  given  in  Reference 
46. 


To  lend  confidence  to  the  analyses  of  the  three  lines  discussed 
above,  in  each  case  an  infinite  array  of  parallel  lines  (G-,  C-,  or 
R-lines),  in  free  space  and  in  the  presence  of  the  flat  earth,  was  an- 
alyzed using  the  same  methodology  but  now  involving  the  transformation 
of  sums  of  cylindrical  waves  into  discrete  spectra  of  plane  waves  (Fourier 
series  expansion).  By  separating  the  lines  in  the  array  sufficiently, 
results  for  the  array  approach  those  of  the  Isolated  line,  providing 
a necessary,  but  not  sufficient,  check  on  the  relevant  mathematics  and 
computer  codes. 

Because  the  formulation  of  the  R-llne  theory  is  long  and  compli- 
cated, rather  than  repeat  it  here,  we  refer  the  reader  to  Reference 
46  and  pass  on  to  a few  numerical  results. 
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B.  Some  Calculated  Results 


For  the  G-line  configuration  shown  in  Figure  55,  curves  of  attenu- 
ation constant  a vs  frequency  are  shown  in  Figures  50  and  57  for  free 
space  and  d=5  meters  over  the  earth,  respectively.  The  phase  constants 
8=w/v  for  the  two  cases  are  Indistinguishable  and  are  shown  in  Figure 
58.  A typical  pattern  of  the  longitudinal  component  of  the  electric 
field,  Ez,  shown  in  Figure  59  on  a polar  plot  with  linear  scale,  dis- 
plays distorting  effect  of  the  ground  at  100  MHz. 


REGION  3Z  ’ 2*0 

<r4  * IO"2 

Figure  55.  Parameters  used  for  numerical  examples. 

The  conductivity  of  copper  was  assumed 
to  be  57  Megasiemens/m. 

For  the  C-line  configuration  shown  in  Figures  60  and  61,  curves 
of  the  attenuation  constant  a vs  frequency  are  shown  in  Figures  62  and 
63  for  free  space  and  d=4  meters  over  the  earth,  respectively.  The 
phase  constants  are  shown  in  Figure  64.  A polar  pattern  of  Ez  is  shown 
in  Figure  65. 
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Figure  57.  Attenuation  constant  a vs.  frequency  of  the  Goubau  line  of  Figure  55 
with  d=5  meters.  The  circled  data  are  calculated  using  the 
theory  of  Chiba  found  in  Reference  51. 


Figure  58.  Phase  constant  8 vs.  frequency  of  the  Goubau  line  of 
Figure  55  with  d=»  (free  space). 
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Figure  59.  E field  distribution  of  Goubau  line  with  d=4  meter  around  a 
concentric  circle  with  radius  r=l  meter.  The  circled  data 
are  those  of  G-array  with  S>10X.  f=100  MHz. 
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Figure  60. 
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Cross-section  of  an  infinitely  long  coaxial  cable 
with  periodic  circumferential  360°  slots 
over  a flat  earth. 
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Figure  62.  Attenuation  constant  vs.  frequency  of  the  C 
of  Figure  60  in  free  space. 
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Figure  63.  Attenuation  constant  a vs.  frequency  of  the  C-llne 
of  Figure  60  with  d«4  meters.  The  circled  data 
are  those  calculated  for  the  C-array 
with  S>!0  . 
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For  the  R-line  configuration  in  free  space  shown  in  Figures  66 
and  67  curves  of  the  attenuation  constant  a vs  frequency  are  shown  in 
Figures  68  and  69  for  the  Goubau-like  (monofilar)  mode  and  the  coaxial- 
like (bifilar)  mode,  respectively.  Figures  70  and  71  show  the  phase 
constants  for  these  two  modes.  These  latter  curves,  if  plotted  on  a 
linear  frequency  scale,  would  be  essentially  linear  with  slopes  cor- 
responding to  v/c  ~ 0.91  for  the  monofilar  mode  and  v/c  ~ 0.75  for  the 
biflar  mode.  Figures  72-75  depict  some  relative  model  field  patterns 
around  the  cable. 

As  In  the  case  of  the  Yagi  line,  it  would  be  most  useful  to  have 
field  plots  normalized  to  one  watt  of  power  transformed  in  a given  mode. 
Although  this  can  be  done,  time  under  this  contract  did  not  permit  it. 

At  the  writing  of  this  report,  data  are  being  generated  to  obtain 
the  mode  propagation  constants  for  an  R-line  near  a lossy  earth  interface. 
Such  data  will  appear  in  Reference  46. 
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Figure  65.  E field  distribution  around  a concentric  circle 
of  radius  r«l  meter  of  a C-llne  of  Figure  60  with 
d=4  meters.  The  circled  data  are  those  calculated 
for  the  C-array  with  S>10A.  f=100  MHz. 
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Figure  66.  Cross-section  of  a coaxial  cable  with  periodic  rectangular 
circumferential  slots  over  a flat  earth. 


o 


X 


Figure  72.  E field  distribution  of  R-line  of  Figure  66  in  free  space 
afound  a concentric  circle  of  radius  r-1  meters. 
f=100  MHz.  Goubau-like  mode. 
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Figure  74.  E field  distribution  of  R-llne  of  Figure  86  In  free  space 
around  a concentric  circle  of  radius  r=l  meter.  f*100  MHz. 
Coaxial-like  mode. 
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Figure  75.  H field  distribution  of  R-llne  of  Figure  66  in  free  space 
afound  a concentric  circle  of  radius  r=l  meter.  f=100  MHz. 
Coaxial-like  mode. 
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IV.  DISCUSSION 

In  this  report  we  have  summarized  theoretical  solutions  for  the 
propagation  constants  and  field  structures  of  surface  waves  supported 
by  two  types  of  infinitely  long  open  transmission  lines  --  the  periodic 
array  of  parallel  rods  or  Yagi  line,  and  the  coaxial  cable  with  periodic 
apertures  or  "leaky  coax".  Both  lines  were  assumed  to  extend  parallel 
to  a lossy  earth  Interface.  In  the  case  of  the  Yagi  line  verifying 
experimental  data  were  obtained  on  a scaled  model  operating  about  2.5 
GHz  and  some  theoretical  data  were  generated  showing  the  transmission 
loss  at  circular  bends  In  the  line  and  reflections  expected  from  a near- 
by (wire  dipole)  intruder. 

Although  a substantial  amount  of  theoretical  data  were  generated 
for  the  Yagi  line,  the  fact  that  so  many  design  parameters  are  avail- 
able (dipole  length,  diameter,  spacing,  height,  frequency,  etc)  pre- 
vented us  in  the  time  alloted  from  calculating  enough  data  to  specify 
a "best"  design.  However,  a few  general  trends  emerged  which  are  useful 
to  anyone  trying  to  use  this  line. 

(1)  In  the  presence  of  a lossy  earth,  the  attenuation  is  a strong 
function  of  dipole  length  (all  other  parameters  being  fixed), 
reaching  a minimum  at  a particular  length.  For  the  range 
of  heights  considered  (1.25X,  1.5X  above  earth)  this  "optimum" 
dipole  length  appears  to  be  about  0.4X-0.44X.  Longer  lengths 
slow  down  the  wave  and  bind  it  closer  to  the  dipole  struc- 
ture, causing  primarily  metallic  losses;  shorter  lengths  speed 
up  the  wave  or  extend  it  so  that  more  attenuation  Is  caused 
by  earth  losses.  It  appears  that  the  wave  should  be  travelling 
about  90*  the  speed  of  light  to  achieve  this  minimum  attenu- 
ation (all  other  things  being  fixed). 
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(2)  All  dipole  lengths  given  in  this  report  are  metric  lengths. 
However,  it  was  found  experimentally,  as  expected,  that  the 
effective  dipole  lengths  are  somewhat  longer  and  they  must 

be  accounted  for  in  accurate  design.  Unfortunately,  the  exact 
increase  in  length  is  not  completely  understood;  it  depends 
upon  rod  diameter,  of  course,  but  it  also  depends  on  the  cur- 
rent distribution  assumed  to  exist  on  the  rods.  In  our  case, 
with  rather  thick  rods  (length-to-diameter  ratio  of  32:1) 
about  a 4%-4%X  increase  in  length  is  appropriate  for  the  piece- 
wise  sinusoidal  current  distribution.  A cosine  on  a pedestal 
current  distribution  demands  an  increase  of  about  3%. 

(3)  The  rod  spacing  also  strongly  controls  speed,  field  extent 
and  attenuation.  We  would  recommend  a spacing  of  about  one 
third  of  a wavelength  in  straight  sections  of  line  and  about 
one  fifth  of  a wavelength  in  curved  sections.  Decreasing 
the  spacing  any  further  on  curves  reduces  shedding  loss  very 
little.  Tapers  in  spacing  at  the  Input  and  output  of  a curve 
were  found  to  effect  negligible  improvement. 

(4)  We  found  that  a reasonably  low  reflection  load  for  the  Yagi 
line  could  be  provided  simply  by  center- loading  several  dipoles 
(we  considered  8 typically)  extending  over  several  wavelengths 
with  resistances  which  linearly  Increase  from  low  to  high 
values  in  the  direction  of  propagation.  We  also  found  that 

a substantial  Improvement  In  radiation  efficiency,  l.e. , a 
reduction  of  radiated  energy  to  the  advantage  of  surface  wave 
energy,  could  be  effected  by  phasing  several  center-fed  dipole 
voltage  sources  In  adjacent  rods  such  that  they  correspond 
to  the  speed  of  the  surface  wave  natural  to  the  structure 
and  induce  constant  currents.  We  also  suggest  the  use  of 
a matched  termination  as  discussed  above  behind  the  phased 
source  to  reduce  radiation  from  the  source  region. 
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(5)  Wire  Intruders  (about  man-high)  encroaching  parallel  to  the 
array  In  free  space  displayed  (for  the  designs  chosen)  reflection 
coefficients  of  the  surface  wave  of  about  0.1  at  an  Intruder 
distance  of  about  1 meter  from  the  line.  This  number  is  nor- 
malized to  the  region  on  the  line  local  to  the  intruder;  back 

at  the  receiver,  the  effect  of  this  reflection  is  much  attenu- 
ated. We  looked  at  the  change  in  impedance  at  the  source 
dipole  due  to  a wire  intruder  and  even  at  very  short  distances 
between  the  two,  the  impedance  change  is  so  slight  as  to  be 
a poor  measure  of  intruder  presence.  Earth  losses  would  de- 
teriorate the  situation  further.  It  appears  that  the  Yagi 
structure  would  not  serve  well  as  a barrier  sensor  of  any 
substantial  length  if  it  is  operated  as  a single  mode  structure. 
Like  most  schemes  which  have  found  some  success  over  long 
distances,  it  may  require  two  parallel  lines,  one  for  trans- 
mitting, one  for  receiving,  to  effect  a large  percentage  change 
in  impedance  upon  introduction  of  an  intruder. 

(6)  Earth  losses  cause  a precipitous  rise  in  attenuation,  becoming 
reasonable  only  for  line  heights  above  the  earth  which  are 
not  practical  at  100  MHz.  A more  practical  version  would 

use  a metallic  ground  screen  which  would  shield  the  wave  at 
least  partially  from  the  earth.  Dipoles  of  length  h=i/2  would 
be  attached  vertically  to  this  screen  and  form  a Yagi  array 
which  could  not  be  penetrated  from  below.  At  100  MHz,  such 
a configuration  would  appear  as  2 foot  high  rods  spaced  about 
a meter  apart. 

(7)  The  Yagi  line  is  quite  dispersive,  so  pulse  distortion  can 
be  expected  if  short  pulses  are  required.  At  100  MHz  a 5% 
bandwidth  corresponds  to  a 200  ft.  pulse,  which  is  an  accept- 
able resolution.  To  minimize  dispersive  effects  over  this 
bandwidth,  it  is  recommended  that  designs  be  used  where  relative 
phase  velocities  exceed  0.80. 
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A theory  of  the  single  coaxial  cable  with  periodic  apertures  was 
completed  during  the  contract  period.  Computer  codes  to  obtain  the 
propagation  constants  and  field  structures  of  the  modes  supported  by 
this  line  were  generated  only  for  the  line  in  free  space;  difficulty 
was  encountered  in  finding  roots  associated  with  propagation  constants 
in  the  presence  of  a lossy  Interface  and  at  the  termination  of  this 
contract  it  was  not  clear  to  us  whether  this  difficulty  was  due  to  poor 
choice  of  starting  values  for  the  iterative  search  scheme  we  employed, 
poor  convergence  behavior  of  the  computational  algorithm,  or  an  error 
in  the  code  itself.  Certain  necessary  conditions,  such  as  consistancy 
of  results  with  and  without  earth  for  large  distances  above  earth,  or 
for  high  frequencies,  were  satisfied.  For  the  free  space  case,  the 
two  modes  on  the  cable  corresponded  to  a Goubau-like  (monofllar)  eigen- 
mode supported  by  the  metallic  shield  and  outer  dielectric  sheath,  and 
a coaxial-like  eigenmode.  The  first  of  these  modes  is  relatively  loss- 
free,  suffering  only  1-2  dB/Km  of  attenuation.  If  however,  the  line 
were  operated  in  the  presence  of  or  buried  in  lossy  earth,  this  mode 
would  be  greatly  attenuated  and  be  rendered  useless  for  intruder  de- 
tection. The  coaxial  mode  suffers  high  attenuation,  of  the  order  of 
150  dB/Km,  even  in  free  space,  but  probably  would  not  suffer  much  more 
in  an  earth  environment.  Most  of  its  field,  however,  resides  inside 
the  coaxial  cable  and  would  be  poorly  coupled  to  an  intruder. 

Future  efforts  should  include: 

(1)  The  study  of  the  Yagi  line  over  a wire  mesh  screen.  The  width 
of  the  screen,  the  mesh  size  and  its  effectiveness  in  screening 
the  lossy  earth  should  be  investigated. 

(2)  The  possibility  of  a multimode  Yagl-structure  should  be  in- 
vestigated so  that  one  mode  could  be  used  for  transmission 
and  the  other  for  reception,  coupling  between  the  two  being 
caused  by  the  asymmetry  introduced  by  an  intruder. 
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(3)  Perhaps  a hybrid  line  could  be  considered,  consisting  of  a 
leaky  coax  attached  along  the  centerline  of  the  Yagl  line. 

Under  symmetry  conditions,  the  Yagl  line  should  not  be  sen- 
sitive to  the  presence  of  the  coaxial  cable  and  so  its  field 
structure  should  be  substantially  unchanged  from  that  studied 
here.  How  severely  an  intruder  (or  imperfections)  would  upset 
this  decoupling  and  cause  a signal  to  be  received  on  the  coaxial 
cable  would  be  a measure  of  the  line's  sensitivity. 

(4)  The  computer  code  to  obtain  the  propagation  constants  of  the 
coaxial  cable  with  periodic  apertures  in  the  presence  of  earth 
should  be  brought  to  completion.  We  simply  did  not  have  the 
time  to  overcome  the  difficulties  encountered  in  this  task 
before  the  contract's  end. 

(5)  The  case  of  two  parallel  leaky  coaxial  cables,  in  free  space 
and  near  earth,  should  be  formulated  since  it  is  this  geometry 
which  has  been  experimentally  successful.  Once  this  theory 

is  developed,  it  makes  sense  to  then  introduce  intruder  ob- 
stacles between  and  above  the  lines  to  estimate  coupling  losses 
and  sensitivity. 

(6)  The  case  of  a leaky  coaxial  cable  along  the  edge  of  a conducting 
septum,  to  represent  such  a cable  running  along  the  top  of 

a chain-link  fence,  should  be  developed.  Since  earth  appears 
so  deleterious,  raising  the  cable  in  this  manner  may  provide 
benefits. 

(7)  Although  not  investigated  at  all  during  this  contract  it  may 
be  profitable  to  consider  a barrier  sensor  which  combines 
the  benefits  of  a relatively  loss-free  transmission  line  with 
the  detection  sensitivity  of  an  antenna.  We  are  thinking 

here  of  a type  of  line  which  guides  the  r.f.  energy  essentially 
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as  a closed  structure,  l.e.,  Insensitive  to  the  environment, 
including  intruders,  from  a source  to  a distant  region  along 
the  line  where  the  line  is  caused  to  radiate,  or  at  least 
extend  its  field.  This  "aperture"  would  be  caused  to  "slide" 
along  the  line  at  a speed  much  slower  than  the  speed  of  r.f. 
propagation,  thereby  performing  a type  of  linear  scanning. 

What  sort  of  structure  is  suitable  to  serve  this  dual  role 
and  how  one  might  signal  the  structure  to  act  as  a closed 
feed  line  or  as  a small  aperture  are  open  questions.  It  cer- 
tainly offers  some  Interesting  possibilities  - for  example, 
an  open-loop  system  could  cause  each  individual  element  in 
a periodic  structure  to  maintain  a radiating  condition  for 
a specified  duration,  at  the  end  of  which  it  would  signal 
its  neighbor  to  do  likewise  much  like  a domino  effect.  A 
more  sophisticated  and  complex  system  incorporating  feedback 
to  the  signalling  mechanism  could  slow  or  stop  the  sliding 
aperture  to  track  an  intruder.  The  concept  is  attractive; 
whether  or  not  it  can  be  made  practical  should  be  investigated. 
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